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I, I ST  OF  SYMBOLS 


h,k,N^,R  Planck,  Boltzman,  Avogadro  and  gas  constants. 

p,T  pressure  and  absolute  temperature 


Q 

* + 
A ,A 


E (= 


EvR^ 


E+,E+ 
v’  r 


* 

k CE  ) 


partition  function 

energized  molecule  and  activated  complex 

total  non-fixed  energy  in  the  active  degrees 
of  freedom  of  a given  energized  molecule. 

* 

vibrational  and  rotational  parts  of  E 

total  non-fixed  energy  in  the  active  degrees  of 
freedom  of  a given  activated  complex 

vibrational  and  rotational  contributions  to  F. 

rate-constant  for  energizing  of  reactant 

rate-constant  for  conversion  of  energized  molecule 
to  products,  Lindeman  theory. 

same  rate-constant  with  energy-dependence,  as  in 
RRKM  theory 


k 


2 


rate  constant  for  de-energizing  of  energized 
moiecule 


*1 


critical  energy  for  reaction 

first  order  rate  constant,  - [l/[  A]  )d[  A] /dt 

statistical  factor  or  reaction  path  degeneracy 

partition  functions  for  adiabatic  rotations  in 
A and  A 


‘h 


* * 

N (E  ) 

X 

Z 


partition  functions  for  active  vibrational  degrees 
of  freedom  in  A 

4* 

number  of  quantum  stages  of  A at  all  energies 
up  to  and  including  E . 

* * 

density  of  quantum  states  of  A at  energy  F.  . 
collisional  energy  transfer  efficiency 
second-order  collision  rate  constant 
entropy  of  activation. 
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I.  INTRODUCTION 


Several  recent  models  have  explained  the  effect  of  pressure  on 
ballistic  properties  as  due  to  an  assumed  change  in  the  microscopic 
processes  that  determine  the  extensive  properties  (1,2).  Each  model 
assumes  that  an  initial  unimolecular  decomposition  reaction  is  the 
rate-determining  step  over  the  entire  pressure  regime  peculiar  to  the 
macroscopic  process.  Both  models  assume  that  external  pressure  causes 
a major  change  in  the  rate  of  this  decomposition  reaction,  thereby 
controlling  the  ballistic  phenomena.  This  report  describes  an  analytical 
study  of  the  effect  of  pressure  on  the  thermal  decomposition  of  a ni- 
trate ester,  n-propyl  nitrate.  The  results  of  the  study  demonstrate 
that  the  pressure-dependence  of  this  reaction  rate  cannot  explain  the 
observed  burning  rate  behavior  of  nitrate  esters. 

The  second  section  of  this  report  describes  the  chemical  reaction 
being  investigated  and  the  methods  used  to  perform  the  analytical  study. 
This  research  involves  the  application  of  the  Rice-Ramsperger- Kassel- 
Marcus  (RRKM)  theory  to  predict  the  rate  constant  for  the  decomposition 
reaction  as  a function  of  pressure.  The  third  section  describes  the 
parameters  used  to  describe  the  reacting  system,  and  provides  details  of 
the  kinetics  calculations.  A final  section  provides  a comparison  of  the 
RRKM  results  and  the  observed  burning-rate  behavior,  and  contains 
comments  on  the  validity  of  the  assumptions  that  underly  the  RRKM 
theory.  Several  appendices  provide  information  about  portions  of  the 
calculation  that  deserve  detailed  explanation. 


II.  UNIMOLECULAR  REACTION  RATE  THEORY 

N-propyl  nitrate  was  chosen  for  study  because  of  good  experimental 
data  on  the  thermal  decomposition  reaction.  The  initial  decomposition 
reaction  is  identified  as  the  scission  of  the  nitro-alkoxy  bond  to 
produce  two  free  radicals,  NO^  and  alkoxy.  This  unimolecular  decom- 
position has  been  studied  for  n-propyl  nitrate  using  the  very  low 
pressure  pyrolysis  (VLPP)  method  (3).  An  RRKM  analysis  of  the  tempera- 
ture dependence  of  the  high  pressure  rate  constant  has  been  performed 
but  no  analysis  of  the  pressure-dependence  of  the  unimolecular  reac- 
tion rate  was  made  (3). 

T. or  'TBir  icu  and  F.A.  Williams,  "Mechanisms  for  the  Steady  Defla- 

gration of  Double-Base  Propellants",  12th  JANNAF  Combustion  Meet- 
ing, CPIA  Publication  Number  273,  1976. 

2.  D. J . Pastine,  M.J.  Kamlett  and  S.J.  Jacobs,  "Volume  and  Pressure 
Dependence  of  Some  Kinetic  Processes  in  Explosives",  Naval  Surface 
Weapons  Center/W0L/TR75-149,  1975. 

3.  G.D.  Mendenhall,  D.M.  Golden,  and  S.W.  Benson,  "The  Very-Low- 
Pressure  Pyrolysis  (VLPP)  of  n-Propyl  Nitrate,  tert-Butyl  Nitrate, 

and  Methyl  Nitrite.  Rate  Constants  for  Some  Alkoxy  Radical  Reactions", 
Int.  J.  Chem.  Kinetics,  7,  725  (1975). 
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The  present  research  consists  of  an  RRKM  analysis  of  the  pressure 
dependence  of  the  nitrate  ester  bond-scission  reaction.  The  follow- 
ing paragraphs  outline  the  development  of  RRKM  theory.  More  complete 
discussions  can  be  found  in  books  by  Forst  (4),  and  Robinson  and 
Holbrook  (5). 

A unimolecular  decomposition  reaction  involves  a single  molecule 
that  undergoes  transformation  to  a reactive  complex  by  some  energization 
process.  The  formation  of  the  complex  is  a function  of  intermolecular 
energy  transfer,  often  via  collision,  and  intramolecular  energy  trans- 
fer after  the  collision  occurs.  Thus,  the  process  is  envisioned  as 

kl  * ka  + 

A s**  A — * A ^products.  (1) 

Most  modern  theories  that  describe  unimolecular  decomposition  reactions 
are  based  on  the  Lindeman  mechanism  (6)  which  assiimes  that  collisions 
between  two  reactant  molecules  will  produce  some  energized  molecules 
with  enough  energy  to  react.  If  conversion  of  reactant  molecules  to 
energized  molecules  is  slow  compared  to  the  rate  at  which  collisional 
de-energization  occurs,  an  equilibrium  population  of  energized  mole- 
cules will  exist  with  concentration  proportional  to  the  concentration 
of  the  reactant  molecules.  Thus  the  rate  constant  will  be  first-order 
with  respect  to  reactant  concentration.  At  low  pressure,  however, 
collisions  cannot  maintain  an  equilibrium  concentration  of  energized 
intermediates,  so  the  rate  constant  will  depend  on  the  rate  of  ener- 
gization which  is  proportional  to  the  square  of  the  concentration  of 
reactant  molecules.  Lindeman' s mechanism  predicts  that  at  low  pressure 
the  rate  constant  increases  with  increasing  pressure,  but  at  high  pres- 
sure the  rate  constant  reaches  a limit  and  remains  constant.  Although 
these  predictions  are  qualitatively  correct,  numerical  studies  based  on 
Lindeman' s mechanism  yield  poor  agreement  between  theory  and  experiment. 

Hinshelwood  extended  Lindeman' s theory  and  laid  the  basis  for  the 
development  of  RRKM  theory  by  suggesting  that  one  needs  to  consider  the 
internal  degrees  of  freedom  when  determining  whether  a molecule  possesses 
sufficient  energy  to  react  (7).  The  chance  for  a molecule  to  have  enough 
energy  to  react  increases  with  the  number  of  degrees  of  freedom  of  the 
reactant  molecule,  primarily  because  of  the  contribution  of  the  vibra- 
tional degrees  of  freedom. 


4.  W.  Forst,  Theory  of  Unimolecular  Reactions,  Academic  Press,  New 
York,  (1973). 

5.  P.J.  Robinson  and  K.A.  Holbrook,  Unimolecular  Reactions,  Wiley- 
Interscience,  New  York  (1972). 

6.  F.A.  Lindeman,  title  not  available,  Trans.  Faraday  Soc. , 1_7, 

S98  (1922). 

7.  C.N.  Hinshelwood,  "The  Theory  of  Unimolecular  Reactions",  Proc . 
Roy. Soc,  (A),  113,  230  (1927). 
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The  final  important  assumption,  suggested  simultaneously  by  Rico 
and  Ramsperger  (8,9)  and  Kassel  (10),  is  that  the  rate  of  decomposition 
of  the  energized  molecule  is  a function  of  the  energy  of  the  species. 
The  species  will  not  decompose  unless  it  possesses  a minimum  internal 
energy,  called  the  critical  energy  Eq,  and  will  react  at  a rate  that 

increases  as  a function  of  the  energy  it  contains  in  excess  of  Eo>  The 
corollary  of  this  is  that  if  energy  is  the  necessary  and  sufficient 
condition  for  reaction,  it  is  implied  that  energy  is  exchanged  among 
the  internal  degrees  of  freedom.  That  is,  there  is  intramolecular 
energy  transfer.  The  developers  of  RRKM  theory  ignore  the  dynamics  of 
such  transfer,  assume  that  the  transfer  is  a random  process,  and  go  on 
to  develop  a statistical  theory  couched  in  the  language  of  transition- 
state  theory  (11). 


The  reaction  scheme  for  the  RRKM  model  is 


A + M 


5k  (E*-*-E*+SE*) 


+ * * * 
A (E  -+E  + oE  ) 


A (E  ) 


k (E  ) 
a 


*■  A 


->  products, 


(2) 


(3) 


* * -k 

The  notation  E ->-£  +5E  means  that  the  energized  molecule  has  energy 

* ★ ★ 

in  the  range  E to  E +6E  . The  second  equation  represents  the  major 
feature  of  the  RRKM  theory  wherein  Absolute  Rate  Theory  is  applied  to 

* 

the  calculation  of  k (E  ) . Thus  we  distinguish  between  an  energized 
* a •* 
molecule  A and  the  activated  complex  A . 


8.  O.K.  Rice  and  H.C.  Ramsperger,  "Theories  of  Unimolecular  Reactions 
at  Low  Pressures",  J.  Amer.  Chem.  Soc.,  4£,  1617  (1927). 

9.  O.K.  Rice  and  H.C.  Ramsperger,  "Theories  of  Unimolecular  Reactions 
at  Low  Pressures  - II",  J.  Amer.  Chem.  Soc.,  5j),  617  (1928). 

10.  L.S.  Kassel,  Kinetics  of  Homogeneous  Gas  Reactions,  Chemical 
Catalog  Co.,  New  York,  1932. 

11.  R.A.  Marcus,  "Unimolecular  Dissociations  and  Free  Radical  Recom- 
bination Reactions",  J.  Chem.  Phys.,  2£,  359  (1952). 
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The  energized  molecule  A is  a molecule  that  contains  in  its 

H * 

active  degrees  of  freedom  a non- fixed  energy  E greater  than  the 
critical  energy  E0.  This  non-fixed  energy  may  consist  of  both 

vibrational  and  rotational  energy. 

The  activated  complex  A+  is  a species  that  is  recognizable  as 
intermediate  between  reactants  and  products,  and  is  characterized  by 
a configuration  that  corresponds  to  a position  atop  a potential  + 
barrier.  The  total  non-fixed  energy  of  the  activated  complex  is  F. 

* 

which  is  equal  to  E less  the  critical  energy  EQ. 


The  RRKM  theory  provides  the  following  mathematical  expression 
for  the  overall  unimolecular  rate-constant,  k 


k 


uni 


VE  )dklfE*+  E%6E*)/k2 


l+ka(E  )/k2(M] 


(4) 


Using  statistical  mechanics,  expressions  for  6kj  an<^ 

k (E  ) are  obtained  leading  to 

cL 


uni 


hQ^ 


oo 

/ 


{£p(E*){  exp(-E*/kT)dE* 


E =E_ 


%-R 

l+kafE*)/ka[  M] 


(5) 


with 


ka(E  ) - 


L*qt 


hQ^N  (E  ) 


E 

£ 

EvR=° 


P<«W* 


where 


C6) 


rotational  partition  function  for  the  activated  complex 
Q,  rotational  partition  function  for  the  reactant 

A 

Q9  vibrational  partition  function  for  the  reactant 
h Planck's  constant 

[P(F.vR)  number  of  quantum  states  of  A at  all  energies  up  to  E 

* * + * 

N (E  ) density  of  quantum  states  of  A at  energy  E 

L*  statistical  factor  or  reaction  path  degeneracy. 


#Non-fixed  energy  is  internal  energy  available  for  redistribution. 


)0 


In  practice  the  product  k-,[M]  is  replaced  by  the  collisional  deactivation 
rate, 

k2[M]  - X Z p (7) 

where 


X,  collisional  efficiency 
Z,  collision  number 
p,  pressure. 


The  RRKM  integral  is  rewritten  as  a sum,  resulting  in  the  follow- 
ing working  equation  for  the  evaluation  of  k 


uni 


LQ\ 

KQjQT 


exp(-HQ/RTJ  AE 


uni 


max 

£ 

L=1 


fIP(E;R)}  exp(-E+/RT) 


l+ka(E  )/XZp 


(8) 


where  the  energy  relations, 

E*  = E + E+ , E+  = (i  ♦ E.  (9) 

o i 

have  been  used.  The  rate  constant  will  have  units  sec"  if  the  ener- 
gies are  expressed  in  kcal  mol  1 and  the  constants  h and  R have  the 
values 


9.S37  * 10 


-14 


kcal  mol"1  sec 


(10) 


R * .0019872  kcal  mol"1  K'1 


(11) 


The  expressions  given  in  equations  6 and  8 are  the  expressions  used 
in  our  study  of  the  decomposition  of  n-propyl  nitrate.  In  the  follow- 
ing section  the  methods  used  to  evaluate  variables  in  these  expressions 
are  reviewed  and  the  model  used  for  the  n-propyl  nitrate  calculations 
is  specified. 


III.  RRKM  PARAMETERS  FOR  N-PROPYL  NITRATE 

There  are  a number  of  variables  that  must  be  determined  before 
a prediction  of  k . can  be  made.  The  computation  of  these  variables 
requires  input  information  about  the  structure  and  properties  of  the 
reactant  molecule  and  the  activated  complex.  In  general,  information 
about  the  reactant  is  available,  whereas  information  about  the  complex 
is  unknown.  Important  data  that  are  required  include: 
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1 

l- 

F 

1.  Vibrational  frequencies  <j 

2.  Moments  of  inertia 

3.  Collision  number  ,3 

4.  Critical  energy,  Eq.  ;j 

Vibrational  frequencies  for  the  reactant  are  needed  for  calcu-  i, 

lation  of  the  vibrational  partition  function,  Q-,  and  the  density  of  H 

* * ^ I 

quantum  states,  N (E  ) . Vibrational  frequencies  for  the  activated  ;j 

complex  are  needed  in  order  to  calculate  the  sum  of  quantum  states  of  J| 

A+.  ji 

l 1 

( s 

The  vibrational  frequencies  for  N-propyl  nitrate  are  known,  but  j; 

those  for  a transition  state  arc  neither  known  nor  measurable.  Often,  i . 

vibrational  frequencies  for  the  transition  state  can  be  assigned  by  j 

transfer  of  data  from  molecules  with  analogous  structures.  In  the  ; 

case  of  bond  stretching  frequencies,  bond  force  constants  arc  predicted  i 

on  the  basis  of  bond-order  correlations  [12).  However,  there  is  no  >, 

well-established  method  for  estimating  the  force  constants  for  bending  ; 

and  torsional  modes. 

Benson  and  O'Neal  have  tabulated  frequencies  for  stretching  and 
bending  modes  and  used  these  values  for  predicting  the  high-pressure  ' 

A-factors  of  unimolecular  reactions  [13,14). 

Any  attempt  to  deduce  the  properties  of  the  activated  complex 
leaves  some  of  the  variables  undetermined.  Generally,  these  variables 
f,  re  used  as  adjustable  parameters  to  make  the  model  predict  the 
correct  value  for  the  experimental  high-pressure  frequency  factor. 

Mendenhall,  Golden  and  Benson  (3)  have  made  no  effort  to  deduce  the 
properties  of  the  complex  since  it  is  known  that  the  calculated  fall- 
off  behavior  is  insensitive  to  the  model  used  provided  that  it  pre- 
dicts the  correct  value  of  A^.  The  high-pressure  rate  constant,  k^, 

is  given  in  RRKM  theory  by  the  expression, 
kT  Q*Q* 

k [RRKM)  = L — -=-=•  exp(-E  /kT)  [12) 

h Q1Q2  0 


12.  H.S.  Johnston,  Gas-Phase  Reaction  Rate  Theory,  Ronald,  New  York, 
1966.  ~~ ~ ~~ 


This  implies  that  the  logarithm  of  the  high-pressure  frequency  factor 
is  given  by 


* *1*2 

An  A (RRKM)  « An  (L  e kT/h)  + d ( T An  ) /dT 

" Wl42 

(33) 

t 

. AS 

= An (L  ekT/h)  * /R  . 

(14) 

A convenient  numerical  value  of  this  expression  is 

l0«1(A.  ■ 10-755  * 1“8l0(L*T)  * iS*/4.756- 

(15) 

In  practice,  the  value  of  AS  is  adjusted  so  that  the  high-pressure 
frequency  factor  agrees  with  the  experimental  value.  The  change  in 
entropy  is  divided  into  vibrational  and  rotational  contributions; 

ASv  = Sv(A+)  - Sv(A) 

(16) 

AS*  - Sr(A+)  - Sr(A)  . 

(17) 

To  solve  for  these  contributions  it  is  assumed  that  the  vibrations  are 
quantum  harmonic  oscillators  and  the  rotations  are  classical  . Thp 
statistical  mechanical  expression  for  the  entropy  is 

„ n „ „ d An  Q 

S = R AnQ  + RT  — . 

(18) 

The  appropriate  expressions  to  calculate  the  vibrational  and  rotational 
entropy  are 

Qv  = [l  - exp(-  h v./kT)]'1  (19) 


An  Qv.  (hv./kT) 

T — i (20) 

[exp(hvi/kT)-l] 

Qr  = 1r’i(8AT/h)3/2(IaIbIc)ls  (21) 

dAn  Q 

T - -d— r-  = 3/2  . (22) 
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Given  these  expressions,  the  use  of  the  RRKM  method  to  study  uni- 
molecular  decomposition  kinetics  is  straightforward.  The  molecular 
parameters  chosen  by  Mendenhall,  et  al,  are  given  in  Table  1 with  one 
exception.  In  this  study,  the  RO-NC^  distance  was  reduced  by  .06  A 

in  both  the  reactant  and  the  complex.  This  change  has  almost  no  effect 
on  the  overall  calculation.  The  parameters  listed  yield  a value  for 
the  logarithm  of  the  high-pressure  A-factor  of  16,5.  This  relatively 
high-value  is  expected  because  the  nitro-moiety  becomes  a free-rotor 
in  the  transition  state. 

In  order  to  complete  the  RRKM  integration  as  a function  of 
pressure,  one  must  specify  the  collision  number,  Z,  and  the  collisional 
efficiency,  X.  The  experimental  apparatus  used  by  Mendenhall,  et  al, 
(3),  is  claimed  to  have  a collisional  efficiency  near  one.  The 
collision  number,  Z,  is  given  by  the  expression 


Z = 

with  Z in 


( a2d  Na/R)  (8tt  Na  k/u^Cl/T)^, 
units  Torr-*  sec-*  when: 


(23) 


0^  = collision  diameter  (cm) 
u = reduced  molar  mass  (g  mol"*) 

T = temperature  (K) 

Na  * 6.0225  * 1023  mol-1 
R = 6.2362  * 104  cm3  Torr  K-*  mol-1 
k a 1.3805  * 10-*6  erg  K-*. 

The  temperature  used  in  this  study  was  600°K,  a temperature  assumed  to 
be  the  surface  temperature  of  a burning  double-base  propellant.  The 

collision  diameter,  o^,  was  chosen  to  be  10-/  cm,  slightly  greater 

than  the  distance  between  the  most  separated  atoms  in  the  equilibrium 
structure  of  n-propyl  nitrate. 
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IV.  RESULTS  OE  THE  RRKM  INTEGRATION  AND 
COMPARISON  WITH  BURNING-RATE  DATA 


The  data  given  in  Table  1 have  boon  used  to  estimate  values  for 

ka (E  ) and  £p(E^)  using  the  Hase-Bunker  computer  code  (IS).  These 

results  and  the  collision  number  have  been  used  to  evaluate  k . as 

uni 

a function  of  pressure  - using  a desk-top  computer  (Tektronix  4051) 
to  calculate  the  rate  constant.  Table  II  lists  the  value  for  as 

a function  of  pressure,  and  the  results  are  shown  in  Figure  1.  It  is 
apparent  from  the  figure  that  the  unimolecular  rate  constant  is  a 
constant  for  pressures  above  ~ 400  Torr,  or  slightly  more  than  one-half 
atmosphere.  Thus,  this  research  predicts  that  the  rate  constant  for 
the  initial  scission  reaction  of  this  medium-size  nitrate  ester  is 
independent  of  external  pressure  for  pressures  above  one  atmosphere, 
for  temperatures  comparable  to  those  at  the  surface  of  a burning  pro- 
pellant. What  does  this  imply  in  regard  to  the  burning  of  nitro- 
cellulose? 

Nitrocellulose,  a polymeric  material,  possesses  an  indeterminately 
large  number  of  vibrational  frequencies.  It  is  a characteristic  of 
unimolecular  reactions  that  as  molecules  become  large,  the  effect  of 
pressure  on  the  rate  constant  disappears;  that  is,  the  fall-off  region 
moves  to  lower  pressure.  Consequently,  we  expect  that  the  effect  of 
pressure  on  k . for  nitrocellulose  decomposition  is  negligible. 

Note,  however , "that  we  mean  only  the  initial  scission  of  the  nitrate 
ester;  our  results  cannot  be  extended  to  the  decomposition  reactions 
of  secondary  radicals  formed  in  the  overall  chemical  kinetic  scheme  of 
nitrocellulose  combustion.  The  next  roport  in  this  series  will 
demonstrate  the  crucial  role  played  by  pressure-variation  for  the  de- 
composition of  small  radicals. 

Given  this  demonstration  of  the  virtual  pressure- independence 
of  the  nitrate-ester  bond-scission  reaction,  it  can  be  concluded  that 
the  pressure-dependence  of  propellant  burning-rate  cannot  be  explaiied 
on  the  basis  of  this  reaction.  This  conclusion,  however,  must  be 
tempered  by  a consideration  of  the  assumptions  made  in  developing 
RRKM  theory,  assumptions  that  limit  the  applicability  of  the  theory. 
There  are  two  assumptions  made  in  the  statistical  development  of  RRKM 
theory. 

The  first  assumption  is  a corollary  of  the  assumption  that  energy 
is  the  necessary  and  sufficient  condition  for  reactions,  it  is 
that  the  internal  energy  of  the  energized  molecule  is  randomly  dis- 
tributed. If  the  energy  is  not  randomly  distributed,  then  the  high 


15.  W.L.  Hase  and  D.  Bunker,  RRKM,  Program  No.  234,  Quantum  Chemistry 
Program  Exchange,  University  of  Indiana,  Bloomington,  Indiana. 
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log  P,  torr 


Figure  1 


Plot  of  the  unimolecular  rate  constant,  k,  . , versus  log  of 
pressure  for  N-propyl  nitrate. 


Table  I.  Molecular  Parameters  for  N-Propyl  Nitrate. 


Vibrational  frequencies, 

Molecule 

Complex 

cm  * , with  degeneracy 

3000 

(7) 

3000  (7) 

in  parentheses 

1600 

1600 

1450 

(4) 

1450  (4) 

1350 

(4) 

1350  (4) 

1200 

1150  (5) 

1000 

1000  (2) 

1150 

(4) 

720  (2) 

750 

(3) 

350 

(4) 

250 

(3) 

125  (6) 

100 

50  (3) 

50 

(2) 

R0-N02,  bond  length, A° 

1.407 

2.407 

2 

Moments  of  Inertia,  amu  A0 

58.1 

68.5 

Xb 

412.3 

542.8 

Ic 

460.8 

601.7 

^red 

39.1 

Table  II.  Rate  Constant  versus  Pressure  Data  for  600°K  Decomposition 
of  N-Propyl  Nitrate. 


Pressure  (torr) 

k . 
uni 

0.5 

24.858 

1. 

30.559 

10. 

44.088 

20. 

45.902 

100. 

47.615 

200. 

47.849 

400. 

47.969 

1000. 

48.04 

10000. 

48.08 

100000. 

48.09 
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pressure  rate  constant  will  not  be  constant,  but  will  show  an  increase 
over  a finite  range  of  pressure  above  the  limit  for  the  decomposition 
of  molecules  with  randomly  distributed  internal  energy.  The  effect 
of  non-random  distribution  of  energy  has  been  discussed  by  Oref  (16)  and 
commented  on  by  Adams  and  Ward  (17).  It  is  unlikely  that  such  an 
effect  will  be  observed  in  molecules  such  as  n-propyl  nitrate  or 
nitrocellulose,  since  the  temperature  remains  relatively  low  during 
the  decomposition  process. 

The  second  assumption  made  in  developing  RRKM  theory  is  that  a 
single  vibrational  potential  defines  the  reaction  hypersurface.  This 
means,  in  dynamic  terms,  that  the  entire  reaction  proceeds  along 
a pathway  corresponding  to  a single  potential  energy  surface.  The 
problem  can  be  understood  with  the  aid  of  Figure  3,  a schematic  that 
shows  the  electronic  energy  for  the  first  three  electronic  states  of 
a nitrate  ester  plotted  as  a function  of  the  alkoxy-nitro  bond  length. 
The  curve  labeled  Wj  corresponds  to  the  energy  of  the  electronic 

ground  state,  while  and  label  the  first  two  triplet  excited 

states.  (This  figure  is  based  on  assignments  first  discussed  by 
L.E.  Harris  (18)).  Curves  Wj  and  W 2 are  everywhere  separated  by  an 

energy  difference,  so  the  state  represented  by  W2  cannot  influence  the 

decomposition  of  a nitrate  ester  in  state  W1 . (It  is  expected  that 

all  the  nitrate  ester  molecules  are  initially  in  the  state  represented 
by  Wj).  If  only  states  and  W2  occur,  then  the  nitrate  esters  meet 

the  second  criterion  for  the  applicability  of  RRKM  theory.  However, 
if  the  state  intersects  the  state  Wj,  as  indicated  in  the  figure, 

then  the  molecule  has  a curve-crossing  in  the  reaction  channel.  In 
this  case,  the  intersection  involves  states  of  different  spins,  and 
coupling  between  such  states  is  said  to  be  "spin  forbidden".  To  the 
cogniscenti  forbidden  is  a hollow  expletive,  since  many  examples  of 
"spin-forbidden"  reactions  proceeding  at  measurable  rates  exist.  The 
curve-crossing  shown,  if  it  is  real,  will  have  the  effect  of  accelerat- 
ing the  rate  of  reaction  at  high-temperature  and  pressure.  In  this 
case,  the  pressure  dependence  of  the  nitrate  ester  decomposition  cannot 
be  determined  by  RRKM  theory. 

The  question  remains,  however,  does  the  curve  crossing  exist? 

Harris  (18)  postulated  the  crossing  on  the  basis  of  the  diffuse  nature 
of  the  electronic  spectra  of  nitrate  esters,  an  interpretation  then 


16.  I.  Oref,  "Nonrandomization  of  Energy  in  Unimolecular  Reaction.  The 
Effects  of  High  Pressures  on  the  Apparent  Rate  Constant",  J.  Chem. 
Phys.,  63,  3138  (1975). 

17.  G.F.  Adams  and  J.R.  Ward,  "Assessments  of  Recent  Models  for  Pressure 
and  Volume  Effects  on  REr.ctions  of  Ballistic  Materials",  Tech.  Rpt. 
ARBRL-TR-0205S,  April  1970.  (AD  #A054640) 

18.  L.E.  Harris,  "The  Lower  Electronic  States  of  Nitrite  and  Nitrate  Ion, 
Nitromethane,  Nitramide,  Nitric  Acid,  and  Nitrate  Esters",  J.  Chem. 
Phys.,  58,  5615  (1973). 
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considered  quite  reasonable.  Recently,  however,  a theoretical  analysis 
by  Pack  (19)  demonstrated  that  there  are  alternative  explanations  of 
diffuse  structure  in  electronic  spectra.  This  being  the  case,  it  is 
not  clear  that  such  a crossing  occurs  in  the  electronic  structure  of 
nitrate  esters.  Experimental  (20)  and  theoretical  (21)  investigations 
of  this  problem  are  underway  in  this  laboratory,  and  a resolution  of 
the  problem  is  anticipated. 

In  any  event,  such  an  esoteric  perturbation  on  the  decomposition 
rate  is  outside  the  ranges  of  validity  of  the  published  models  for 
the  effects  of  pressure  on  ballistic  properties,  since  both  models 
were  limited  to  consideration  of  the  effects  of  pressure  on  vibra- 
tional frequencies.  It  is  concluded,  on  the  basis  of  the  analysis 
described  in  this  report,  that  the  initial  decomposition  reaction 
of  the  nitrate  esters  contained  in  double-base  propellants  cannot  be 
the  rate-determining  step  for  the  observed  variation  of  burning-rate 
with  pressure. 


19.  R.T.  Pack,  "Simple  Theory  of  Diffuse  Vibrational  Structure  in 
Continuous  UV  Spectra  of  Polyatomic  Molecules.  I.",  J.  Chem.  Phys., 
65,  4765  (1976). 

20.  R.  Beyer,  private  communication. 

21.  G.F.  Adams,  G.D.  Purvis  and  R.J.  Bartlett,  unpublished  calculations 
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APPENDIX  A.  COMMENTS  ON  DETAILS  OF  THE  RRKM  CALCULATION 


The  working  expressions  for  k ^ and  kft  (E  ) are  repeated  here  for 


convenience. 


lV*  "ax 

ku„i  ‘ hQ^  “P<-VRT)  4E*  J 


(£P(E*R)}  exp(-E+/RT) 


1 + k (E  )/AZP 

a 


(A.  1) 


ka(E  ) 


l\+  V 

hQiN  (E  } h;r=o 


(A.  2) 


For  the  sake  of  completeness,  the  methods  used  to  evaluate  the  sums  and 
densities  of  molecular  quantum  states  should  be  outlined.  Application 
of  the  RRKM  method  is  critically  dependent  in  the  methods  used  to 
evaluate  the  density  of’quantum  states  N*(E*R)  for  the  reactant  A and 

the  sum  of  quantum  states  ZP(E+„)  for  the  postulated  activated  complex 
+ * VK  # 

A . The  computation  of  N and  £P  for  n-propyl  nitrate  was  performed 
using  the  semiclassical  method  developed  by  Whitten  and  Rabinovitch  (1), 
and  this  method  will  be  described  briefly. 

The  density  of  quantum  states  N(E)  can  be  obtained  from  the  sum  of 
quantum  states  W(E)  = EP(E)  by  differentiation  with  respect  to  energy. 
That  is 

N(E)  ,d»|Ei=diiPiaij  (A.S) 

so  we  need  only  find  a differentiable  form  for  W(E).  For  a set  of  s 
oscillators,  each  set  of  quantum  numbers  v^  gives  rise  to  one  quantum 

state.  The  total  number  of  states  at  nil  energies  up  to  and  including 
Ev  is  given  by, 


# We  shall  hereafter  adapt  the  notation  W(E)  = £ P(E),  and  omit 
superscripts . 

1.  G.Z.  Whitten  and  B.S.  Rabinovitch,  "Accurate  and  Facile  Approxi- 
mation for  Vibrational  Energy-Level  Sums",  J.Chem.  Phys.  38,  2466 
(1963) . 
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(A. 4) 


l pey 

E *0 
v 


l 1 •••!  i 


s=0 


Vr° 


V° 


where  the  summation  is  taken  over  all  such  that  the  total  vibrational 
energy  is  less  than  or  equal  to  E . Inthe  classical  approximation, 
this  sum  is  replaced  by  an  integral  over  the  occupation  number,  v^. 

The  evaluation  of  the  integral  is  given  by  Robinson  and  Holbrook  1f2), 
in  Appendix  6,  and  the  result  is 


W(EV)  - P(Ey)  * 


s.'nhv. 


(A.  6) 


Marcus  and  Rice  (3)  extended  the  classical  treatment  by  suggesting  that 
an  oscillator  with  non-fixed  energy  has  a total  vibrational  energy 

above  the  ground  state  of  Ev+Ez,  where  is  the  zero-point  vibrational 

energy.  Modification  of  the  classical  expressions  (A. 5)  and  (2a)  yields 

, s 

(A. 7) 


W(EV) 


(E  +E  )' 
v v z 

s!  Hhv. 

1 


NCEJ  = 


(E  +E.) 


s “ 1 


(A. 8) 


nhv. 


These  expressions,  though  more  accurate  than  the  classical  approximations, 
always  overestimate  the  number  of  states  at  any  level . An  improved 
version  of  the  semiclassical  expressions  has  been  developed  by  Rabinovitch 
and  co-workers. 


It  was  suggested  that  only  a fraction  of  the  zero-point  energy 
should  be  added  to  the  vibrational  energy  E . An  empirical  factor  a 
was  introduced  to  multiply  E2  in  equations  V(A.7)  and  (A. 8).  Whitten 
and  Rabinovitch  discovered  a relation  between  the  parameter,  a,  and  a 
modified  frequency  dispersion  parameter,  8,  defined  by 


6 


s-1 

s 


(A.  9) 


The  relation  is 

/ 

a - 1 - 6 • u(E  ) (A. 10) 


2.  P.J.  Robinson  and  K.A.  Holbrook,  Unimolecular  Reactions,  Wiley- 

Interscience,  New  York  (1972).  ” 

3.  R.A.  Marcus  and  O.K.  Rice,  "The  Kinetics  of  the  Recombination  of 
Methyl  Radicals  and  Iodine  Atoms",  J.  Phys.  and  Colloid.  Chem. , 
55,  894  (19S1) . 
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where  E; 


and  w(E/)  is  a function  of  E'  defined  by 


(0.1  < E;  < 1.0) 


a)  = (S.  E;  + 2.73  E/0,S  + 3.51)'1 2 3 4 5 6 


(A. 11a) 


(1.0  < E'  < 8.0) 


u>  = exp(-2.4191E/  0,2S) 


(A. lib) 


The  molecular  frequencies  are  used  to  calculate  $ and  E^,  and  then  for 

a non-fixed  energy  Ev>  the  reduced  energy  E*  is  calculated  and  used 

to  obtain  ^(E7 8 9 10).  Given  this  data,  the  parameter  a and  the  sum  of  states 
can  be  computed.  The  density  of  states  can  then  be  computed  using  the 
expression 


N(EV) 


(Ev*aE2)s-‘ 
(s-i) inhv. 


where  f is  obtained  by  differentiating  (A. 11a)  and  (A. lib). 

The  Whitten  Rabinovitch  method  has  been  thoroughly  tested,  and  the 
results  are  in  good  agreement  with  the  direct-count  for  non-fixed 
energies  Ey  greater  than  E , and  often  for  lower  energies.  Whitten 
and  Rabinovitch  have  extended  the  treatment  to  include  the  case  where 
active  rotations  need  to  be  taken  into  account.  Again,  the  equations 
suggested  have  been  widely  tested  and  the  results  appear  to  be  good 
approximations  to  W(E  ) and  N(E  ). 


APPENDIX  B:  DATA  FORMAT  FOR  RRKM 

The  computer  code  RRKM  is  stored  as  an  UPDATE  file  on  the  BRL 
Cyber  173.  This  code  can  be  used  by  submitting  the  following  deck  of 
cards , 


1 . JOBCARD 

2.  ACCOUNTING  CARD 

3.  ATTACH  (0LDPL, RRKM, ID-GEOFADAMS) 

4.  UPDATE (F) 

5 . FTN , I . 

6.  LGO. 

7.  *E0R 

8.  * ID, (MARY)  +-  arbitrary  choice  of  ID 

9.  *E0R 
DATA 

10.  *E0I . 


★ 

This  job  will  produce  values  for  N(E  ),  w(E 


and  k (E  ) . 
a 


These  data 
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6 v. 


can  be  used  to  calculate  k ^ as  a function  of  pressure  by  solving 

equation  (A.l).  Unless  one  has  substantial  experience  with  such 
computations,  it  is  recommended  that  the  computation  of  k ^ be  carried 

out  on  a desk-top  computer  or  as  an  interactive  exercise  on  a terminal. 
By  following  this  procedure,  one  need  not  compute  kun^  for  an  excessive 

pressure  range.  Users  are,  however,  invited  to  disagree  with  this 
philosophy . 

The  data  required  to  perform  a standard  RRKM  calculation  is  as 
follows : 

1.  Format  3F  10.0,  IS,  F10.0.  A card  with  Eq  (threshold  energy), 

* 4- 

E (lowest  molecular  energy  used  in  the  calculate  .is),  AE  (upward 
increment  of  same),  n,  the  number  of  successive  molecular  energies  for 
which  the  calculation  will  be  made,  and  a , the  statistical  factor. 

2.  Two  title  cards.  Say  anything  you  like. 

3.  Format  2110,  n^,  the  number  of  reactant  frequencies  and  nr> 
the  number  of  free  internal  rotors. 

4.  Format  7F10.0.  The  molecular  frequencies,  seven  per  card. 

-1  1 
Units,  cm  . 


5.  Reduced  moments  of  inertia  for  the  internal  rotors  of  the 

2 

reactant.  Units  amu  A°  . Format  7F10.0.  Omit  this  card  if  N * 0. 

r 

6.  Symmetry  numbers  corresponding  to  item  S.  Same  format. 

Omit  if  n =0. 

o 


7. 

Format, 

110. 

KROTS0. 

8. 

FORMAT, 

3F.10. 

Three  overall 

moments  of 

inertia  for  the 

reactant 

; units  amu-A*". 

9. 

FORMAT, 

110. 

NANH=0. 

10. 

FORMAT, 

2110. 

JDEN=0,  JSUM=0. 

11. 

FORMAT, 

110. 

NCRIT=0. 

12. 

FORMAT, 

2110. 

♦ , + 

and  n , same  as  item 

3 but  for  the 

critical 

configuration. 

+ + 

n£  + n = N.  + 
f r f 

Nr  - 1. 

13.  Frequencies,  like  item  4 except  for  critical  configuration. 


26 


14.  Reduced  moments,  like  item  5,  except  for  critical  con- 
figuration. 

15.  Symmetry  numbers,  like  item  6,  except  for  critical  con- 
figuration. 

16.  The  overall  moments  of  inertia  for  the  critical  configuration 

as  in  item  8.  1^,  1^,  and  sequence  must  be  the  same  as  for  the 

molecule. 


17.  Format  F10.0.  Rotational  temperature  for  adiabatic 
rotations,  °K. 

18.  Format  110.  NTRAN  = 0. 

Input  may  be  stacked.  To  stop  the  program  add  a blank  card  after 
item  18. 
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